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1  INTRODUCTION 


Many  materials  have  different  elastic  behavior  in  tension  and  compression. 

A  few  examples  of  such  materials  are  concrete,  rock,  tire-cord  rubber,  and  soft 
biological  tissues.  As  early  as  1864,  St.  Venant  [1]  recognized  this  behavior  by 
analyzing  the  pure  bending  behavior  of  a  beam  having  different  stress-strain 
curves  in  tension  and  compression.  Timoshenko  [2]  originated  the  concept  of  bi¬ 
modulus  (or  bimodular)  materials  in  1941  by  considering  the  flexural  stresses  in 
such  a  material  undergoing  pure  bending.  Ambartsumyan  [3]  in  1965  renewed  interest 
in  the  analysis  of  bimodular  materials,  i.e.,  materials  having  different  moduli  in 
tension  and  compression.  Since  then,  there  have  been  numerous  investigations  on 
the  static  behavior  of  bimodular  beams;  these  were  surveyed  by  Tran  and  Bert  [4]. 
Recently,  Bert  and  Gordaninejad  [5]  studied  bending  of  thick  beams  of  "multimodular" 
materials. 

Only  a  few  studies  have  been  made  on  vibration  of  bimodular  beams.  Recently, 
Bert  and  Tran  [6]  worked  on  transient  response  of  thick  beams  of  bimodular  materials. 

The  present  paper  deals  with  the  forced  vibration  of  beams  made  of  "multi¬ 
modular"  materials.  The  transfer-matrix  method  [7],  which  computationally  is  very 
efficient,  is  applied.  Also,  the  beam  is  modeled  as  a  Timoshenko  beam,  i.e.,  both 
transverse  shear  deformation  and  rotatory  inertia  are  considered. 

2  MODELING  OF  THE  STRESS-STRAIN  CURVE 

The  nonlinearity  of  the  stress-strain  curve  is  one  of  the  main  difficulties 
arising  in  structures  undergoing  even  moderate  deflections.  Piecewise  lineariza¬ 
tion  of  the  stress-strain  relation  has  been  applied  to  overcome  this  problem. 

Durban  and  Baruch  [8]  used  a  floating  piecewise  linear  approximation  to  construct 
the  two  "best"  straight  lines  approximating  the  Ramberg-Osgood  stress-strain  relation 
[9].  Bert  and  Kumar  [10]  recently  presented  experimental  stress-strain  curves 
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for  unidirectional  cord-rubber  materials  and  expressed  the  curves  in  tudwik 
power-law  form  with  different  coefficients  and  exponents  in  tension  and  compres¬ 
sion. 

In  the  present  work,  a  stress-strain  curve  for  aramid-rubber  taken  from 
[10]  has  been  linearly  approximated  by  four  segments  (two  segments  in  tension 
and  two  segments  in  compression;  see  Fig.  1).  For  choosing  the  "best"  two  straight 
lines,  the  area  between  two  fitting  lines  and  the  experimental  curve  in  each  por¬ 
tion  has  been  minimized.  To  find  comparable  moduli  for  the  bimodular  case,  one 
has  to  minimize  the  area  between  two  straight  lines  and  the  experimental  curve. 
Finally,  for  the  "unimodular"  case,  the  "best"  single  straight  line  is  used  (see 
Appendix  A) . 


3  CLOSED- FORM  SOLUTION 

Consider  a  solid  rectangular-cross-section  beam  of  thickness  h  and  length 
i.  The  beam  coordinates  are  taken  such  that  the  xy-plane  coincides  with  the  mid- 
plane  of  the  beam  and  the  z-axis  is  measured  positive  downward.  For  a  four- 
segment  approximation  of  the  normal  stress-strain  curve,  considering  the  general 
case  (i.e.,  when  -  h/2  <  ac  ,  at  <  h/2),  the  following  stress  field  has  been 
considered  for  the  case  of  convex  bending  (see  Figs.  1  and  2). 


E1C  elC  +  E2C(£x-£lC) 

ElCcx 

Eltcx 

ElV  +  E2t(ex-elt} 


-  h/2  <  z  <  ac 

ac  <  2  <  2n 

2n  <  2  <  at 
at  <  z  <  h/2 


(1) 


(2) 


where  E^.EgSE^.I^.G.ei6,  and  are  material  constants,  is  the  axial  nor¬ 
mal  stress,  ex  is  the  axial  normal  strain,  yx2  is  the  transverse  shear  strain,  tx2 

Is  the  transverse  shear  stress,  and  z„  is  the  location  of  the  neutral  surface. 

n 


Fig.  1  Multimodular  model . 


wx 


Fig.  2  Stress  Distribution  of  a  Multimodular  Beam  for 
Convex  Bending 
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It  is  noted  that  this  material  is  linear  elastic  in  shear.  Comparison  of  Figs. 
1  and  2  leads  to 


‘x  =  ‘(2-zn> 


'1 


c 


t 

£1 


) 

) 


Using  linear  strain  measure,  one  obtains 


(3) 

(4) 

(5) 


e  =  U  =  U  +  Z\b 

X  ,x  ,x  *,x 

Y  =  W  +U  =  W  +  li; 
rxz  ,x  ,z  ,x  v 


Comparison  of  equations  (3)  and  (6)  gives 


u  =  -*zn 
•  x  n 


Note  that  (  )  denotes  3(  )/3x. 


(6) 


(7) 


Timoshenko  beam  theory  is  implemented  here,  by  using  the  definitions  of  the 
normal  and  transverse  shear  stress  resultants  and  moment,  each  per  unit  width  as 
follows : 


h/2 


h/2 


(N,Q)  = 


(VTxz)dz 


M  = 


zax  dz 


(8) 


-h/2  -h/2 

One  can  write  the  constitutive  relation  for  a  multimodular  beam  as 


V 

rA+ cA 

8tCNB 

0 

r  > 

U,x 

A' 

B' 

0 

r  \ 

u,x 

« 

►  * 

b+4 

D*C» 

0 

'.X 

►  a 

B" 

D' 

0 

4 

♦.X 

0 

0 

s 

w 

X 

+ 

-e- 

0 

0 

S 

w  +  * 

;  ,x  J 

(9) 


where  A,  B,  D,  and  S  denote  the  respective  extensional,  flexural -extensional 
coupling,  flexural,  and  transverse  shear  stiffnesses  defined  by 
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(A,B,D)  *  (l,z,z2)E^k)dz  k=t!c  S  =  K2  G  dz 


ABB  D 

Here,  the  stiffnesses  C^,  CN,  CM,  and  CM  are  not  present  in  unimodular  or  bi~ 

modular  materials  (see  Appendix  B).  In  equation  (10),  t  and  c  denote  tensile- 

2 

strain  and  compressive-strain  regions,  respectively,  and  K  is  the  shear  correction 
coefficient*.  The  general  equations  of  motion,  if  zn  (neutral -surface  location) 
is  constant  along  the  beam  are 


where 


A\xx  * 

B'* 

,xx 

Pu,tt  X-  R»>tt 

S(\xx  * 

H>  J  ■  Pw 

<B"Vx  1 

'  °\xx> 

-  S(w  +  \p)  : 

9  * 

h/2 

f 

(P.R.I)  ■ 

P(1 

,z,z2)dz 

-h/2 

and  p  is  the  density  of  material. 

« 

For  guided-guided  boundary  condition,  i.e., 

u(0,t)  *  u(*,t)  =  0  ;  'I'CO.t)  =  i|>U,t)  =  0 

(12) 

Q(O.t)  =  Q(t.t)  =  0 
if 

q(x,t)  *  qQ  cos  ax  cos  fit  (13) 

then  the  following  sets  of  functions  satisfy  the  equations  of  motion 
u(x,t)  =  U  sin  ax  cos  nt  ;  <pCx,t)  =  ?  sin  ax  cos  nt 

(14) 

w(x,t)  »  W  cos  ax  cos  nt 

*In  actuality,  enforcement  of  the  axial-free  equilibrium  equation  of  elasticity 
requires  that  K*  for  a  multlmodular  beam  (even  a  single-layer  one,  such  as  treated 
here)  be  a  function  of  the  level  of  normal  strain  (through  the  piecewise  segmenta¬ 
tion  of  the  stress-strain  curve).  However,  in  this  paper,  Kz  is  assumed  to  be  a 
constant. 


where 


Q  =  2nf  ,  a  -  mir/z  (m  =  1 ,2,3, . . . )  (15) 

f  =  circular  frequency  of  the  excitation,  and 


q  (Sa)(B'a2-Rfl2) 


(Sa2-Pfi2)[B'a2-Rn2) (B"a2-Rn2)-(Oa2+S-In2)(A'a2-Pn2)]+(Sa)2(A'a2-Pn2) 

(16) 

:  -  A* a2  -  Pfl2  U  ;  u - -  [q  -  (Sa)(A‘a2-Pn2)-j  (j 

B<a2.Rn2  Sa2-Pn2  °  B'a2-Rfi2 


Since  from  equations  (7) 

then 


-  u 


,x 


"n 


=  B'a2-Rfl2 
A’a2-Pn2 


=  constant 


(17) 

08) 


4  TRANSFER-MATRIX  SOLUTION 

The  transfer-matrix  model  used  in  the  present  study  is  the  same  as  that 
employed  in  [5]  except  for  the  station  matrix  (see  Appendix  C),  which  here  in¬ 
cludes  more  terms  due  to  the  motion.  The  transfer  matrix  for  the  assumed  beam 
is  of  the  following  form 

N  •  • 

*-s-*ns+i  =  ij1  awVi^fWz^o  (19) 

where  [T-].  is  the  field  matrix,  [Te].  is  the  station  matrix,  Nc  is  the  number 

T  1  5  l  s 

of  stations,  aji/2  is  the  length  of  each  of  the  half  fields  at  the  ends  of  the 
beam,  m  is  the  length  of  each  of  the  whole  fields,  and  [S]N  +1 ,  [S]Q  are  state 

T  ^ 

vectors,  i.e.,  (u,w,T,N,Q,M)  ,  at  the  two  ends  of  the  beam. 

In  the  calculation  of  the  stiffnesses  for  the  cases  where  the  axial  force  is 

not  zero,  the  neutral -surface  location  and  the  corresponding  distances  to  the 

"break  points"  (ac  and  at)  in  the  ax  vs  z  curve  are  not  constant  and  not  known 

a-priori.  Therefore,  an  iterative  technique  has  been  employed  to  compute  the 

neutral -surface  locations  z„,  also  a„  and  a*.  One  must  first  assume  (2N  +2) 

net  s 
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sets  of  values  of  zn,  ac,  and  at  and  then  compute  the  sitffnesses  and  solve  the 
governing  equations  for  the  state  vector.  Finally,  by  using  equations  (C.l), 

(C.3),  and  (C.4),  compute  the  new  values  of  zn>  ac>  and  at-  Obviously,  if  the 
assumed  and  computed  sets  of  zn,  ac>  and  at  are  in  sufficiently  close  agreement, 
the  problem  is  solved;  otherwise,  assume  the  calculated  set  zn>  ac,  and  at  and 
repeat  the  procedure. 

5  NUMERICAL  RESULTS 

The  numerical  results  are  presented  for  a  thick,  multimodular  beam  with  a 
rectangular  cross  section.  The  material  of  the  beam  is  chosen  to  be  aramid  cord- 
rubber  which  is  used  in  automobile  tires  (see  Table  1).  Four  different  boundary 
conditions  are  investigated  (see  Table  2)  and  comparisons  are  made  between  multi¬ 
modular,  bimodular,  and  unimodular  models  for  each  set  of  boundary  conditions. 

In  this  study,  a  mesh  of  twenty-five  elements  is  used  with  each  element  being  of 
length  of  0.32  in.  The  shear  correction  coefficient  is  taken  to  be  5/6. 

In  order  to  validate  the  transfer-matrix  solution  (TMS),  Fig.  3,  a  compari¬ 
son  is  made  between  the  closed-form  solution  (CFS)  and  the  TMS  for  a  guided-guided 
beam  with  cosine  load  distribution  (case  1).  Also,  a  comparison  is  made  among 
unimodular,  bimodular,  and  multimodular  (static  and  dynamic)  cases  (see  also 
Table  3)  for  f  =  100  Hz.  As  one  can  see,  there  is  excellent  agreement  between 
the  TMS  and  CFS  results.  However,  this  agreement  can  be  improved  even  further 
by  increasing  the  number  of  elements. 

For  the  other  cases  (2-4),  the  CFS  is  not  available;  therefore,  in  Figs.  4, 

6,  and  8  comparisons  between  different  models  (one,  two,  and  four  segment  approxi¬ 
mations)  are  made.  As  one  might  notice  in  all  four  cases,  there  is  considerable 
difference  between  transverse  deflection  of  multimodular  and  bimodular  beams  on 
one  hand  and  that  of  the  unimodular  model  on  the  other  hand.  In  contrast,  there 
is  no  substantial  difference  between  multimodular  and  bimodular  results. 

Another  interesting  observation  in  Fig.  4  is  that  for  f  =  100  Hz,  the 
unimodular  beam  is  in  the  range  of  its  first  mode,  whereas  the  bimodular  and  multi- 


astic  Properties 
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Table  1  Elastic  properties  and  geometric  parameters  for  an 
aramid-cord  rubber  beam 


Longitudinal  Young's  Modulus, 
MPa  (psi  xlO-6) 

Model* 

Tension 

Compression 

Tension  and  Compression 

E,t 

E,c _ EoC 

G 

M 

4000  2896 

(0.580)  (0.420) 

■ 

3.70  (0.537) 

h' 

Ebc 

B 

1 

WHihoW 

124 

(0.018) 

3.70  (0.537) 

E 

E 

U 

■ 

mmmm 

* 

3.70  (0.537) 

Beam  length  20.32  cm  (8.0  in.) 

Beam  thickness  1.52  cm  (0.6  in.) 

Beam  width  2.54  cm  (1.0  in.) 

O  S- 

•jr  <u 
s- 


£  <o 

o 

a; 

13  Q- 


M  'v  multimodular,  B  ^  bimodular,  U  *  unimodular. 


Table  2  Summary  of  cases  considered 


DEFLECTION 


*1 

* 

BIMODULAR  CFS  &  TMS 

o 

UNIMODULAR  CFS 

_  __  _ 

UNIMODULAR  TMS 

j  -  MULTIMODULAR  CFS  &  TMS 

-0.05- 


-0.  10- 


-0.20- 


DIMENSIONLESS  POSITION,  x/t  (FOR  HALF  OF  THE  BEAM) 

Fig.  3  Comparison  among  multimodular,  bimodular,  and  unimodular  deflection 
distribution  for  closed-form  and  transfer-matrix  solutions  of  guided-quided 
aramid-cord  rubber  beam  (f  *  100  Hz) 


Table  3  Comparison  between  CFS  and  TMS  for  an  aramid-cord  rubber  beam 
(case  1 )f,  f  =  100  Hz 


X/H 

N  x  102 

.  lb 

Q  x 

10,  lb 

M, 

1 b-in . 

CFS 

TMS 

CFS 

TMS 

CFS 

TMS 

0.00 

-0.169 

-0.164 

0.000 

0.000 

0.203 

0.203 

0.02 

-0.168 

-0.162 

-0.198 

-0.200 

0.201 

0.202 

0.06 

-0.157 

-0.152 

-0.585 

-0.587 

0.189 

0.189 

0.10 

-0.137 

-0.132 

-0.935 

-0.937 

0.164 

0.165 

0.14 

-0.108 

-0.104 

-1.226 

-1.228 

0.129 

0.130 

0.18 

-0.072 

-0.069 

-1.141 

-1.143 

0.086 

0.087 

0.22 

-0.032 

-0.031 

-1.565 

-1.566 

0.038 

0.038 

0.26 

0.011 

0.010 

-1.590 

-1.591 

-0.013 

-0.013 

0.30 

0.052 

0.050 

-1.514 

-1.516 

-0.063 

-0.063 

0.34 

0.091 

0.087 

-1.344 

-1.346 

-0.109 

-0.109 

0.38 

0.123 

0.119 

-1.089 

-1.091 

-0.148 

-0.148 

0.42 

0.148 

0.143 

-0.766 

-0.768 

-0.178 

-0.178 

0.46 

0.164 

0.158 

-0.394 

-0.396 

-0.196 

-0.197 

0.50 

0.169 

0.163 

0.000 

0.000 

-0.203 

-0.204 

0.54 

0.164 

0.158 

0.394 

0.396 

-0.196 

-0.197 

0.58 

0.148 

0.143 

0.766 

0.768 

-0.178 

-0.178 

0.62 

0.123 

0.119 

1.089 

1.091 

-0.148 

-0.148 

0.66 

0.091 

0.087 

1.344 

1.346 

-0.109 

-0.109 

0.70 

0.052 

0.050 

1.514 

1.516 

-0.063 

-0.063 

0.74 

0.011 

0.010 

1.590 

1.591 

-0.013 

-0.013 

0.78 

-0.032 

-0.031 

1.565 

1.566 

0.038 

0.038 

0.82 

-0.072 

-0.069 

1.141 

1.143 

0.086 

0.087 

0.86 

-0.108 

-0.104 

1.226 

1.228 

0.129 

0.130 

0.90 

-0.137 

-0.132 

0.935 

0.937 

0.164 

0.165 

0.94 

-0.157 

-0.152 

0.585 

0.587 

0.189 

0.189 

0.98 

-0.168 

-0.162 

0.198 

0.200 

0.201 

0.202 

1 .00 

-0.169 

-0.164 

0.000 

0.000 

0.203 

0.203 

*CFS  ^ 

closed-form 

sol ution; 

TMS  -v  transfer-matrix 

solution 

f  +  0.3305  0  <  x/i  <  0.22 

+For  case  1,  Z  =  z  /h  is  piecewise  constant,  equal  to  <  -  0.3305  0.22  <  x/e.  <  0.78 
n  l  +  0.3305  0.78  <  x/l  <  1.00 


DIMENSIONLESS  POSITION,  x/l  (FOR  HALF  OF  THE  BEAM) 

Fig.  6  Comparison  among  multlmodular,  bimodular,  and  unimodular  deflection 
distribution  for  transfer-matrix  solution  of  clamped-clamped  aramid-cord 
rubber  beam  (f  =  100  Hz) 


DEFLECTION 
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Fig.  8  Comparison  among  multimodular,  bimodular,  and  unimodular  deflection 
distribution  for  transfer-matrix  solution  of  clamped-free  aramid-cord  rubber 
beam  (f  ■  100  Hz) 


modular  beams  are  between  their  first  and  second  modes.  The  explanation  for  this 
is  that,  for  this  case  (2),  most  of  the  layers  of  the  beam  are  under  compression 
and  since  E^c,  E2C,  and  E^c  are  much  smaller  than  E  (see  Table  1),  then  the  uni- 
modular  beam  is  stiffer  than  the  other  two.  Therefore,  the  fundamental  frequency 
of  the  unimodular  beam  is  higher  than  those  of  the  bimodular  and  multimodular 
beams. 

Also,  the  distributions  of  u  (axial  displacement),  J  (bending  slope),  N  (axial 
force),  Q  (transverse  shear  force),  and  M  (bending  moment)  are  shown  graphically 
in  Figures  5,  7,  and  9  for  f  taken  to  be  100  Hz.  Note  that  for  cases  1  and  3 
this  frequency  is  less  than  the  fundamental  frequency,  whereas  for  cases  2  and  4, 
the  frequency  of  100  Hz  is  in  the  range  of  the  first  and  second  modes.  The  first 
three  mode  shapes  of  a  clamped-free  beam  of  multimodular  material  is  investigated 
in  Fig.  10.  For  this  case,  the  natural  frequencies  associated  with  the  first 
three  modes  are  f-j  *  30.9  Hz,  f2  *  131.1  Hz,  and  f3  *  278.5  Hz. 

Finally,  by  rewriting  the  equations  of  motion  in  a  new  form  (N  x  =  P1  u  tt, 

M  -  Q  s  I*  P  ’  -  q(x,t)),  the  effect  of  translatory  and  rotatory 

inertia  coefficients  on  axial  force  for  a  thick  multimodular  clamped-clamped 
beam  (f  *  100  Hz)  is  studied  (see  Table  4).  The  results  show  the  significant 
effect  of  I  and  slight  effect  of  as  one  looks  at  it  through  full  theory 
(vibration)  as  compared  to  the  static  case  (I  *  P.|  =  P2  a  0) . 

6  CONCLUSIONS 

An  analysis  of  forced  vibration  of  a  thick  beam  with  a  rectangular  cross 
section  and  made  of  "multimodular"  material  is  presented.  In  this  study, 
numerical  results  obtained  by  both  the  closed-form  and  transfer-matrix  methods 
are  given  for  a  beam  made  of  aramld-cord  rubber. 

Comparisons  are  made  on  one  hand  between  closed-form  and  transfer-matrix 
results  and  on  the  other  hand,  among  unimodular,  bimodular,  and  multimodular 


MAXIMUM  DEFLECTION,  wU)  x  10 


FREQUENCY,  f,  Hz 


Fig.  10  Investigation  of  the  first  three  mode  shapes  of  thick  multimodular, 
clamped- free  aramld-cord  rubber  beam  with  rectangular  cross  section. 
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Table  4  Effect  of  translatory  and  rotatory  inertia  coefficients  on  axial 
force  for  a  thick  multimodular  cantilever  beam  (f  =  100  Hz) 


Axial  force,  lb  x  103* 

x/L 

I.Pl.Pztf) 

Full  Theory 

1=0  Pi&P2*0 

I&P !=0  P2f0 

P ! =0  I&P2/0 

I,Pi,P2=0 

Static 

11.74 

11.76 

-0.1399 

-0.0699 

11.49 

11.51 

-0.1399 

-0.0699 

10.56 

10.58 

-0.1399 

-0.1621 

-0.0699 

0.10 

8.90 

8.92 

-0.1399 

-0.1621 

-0.0699 

0.14 

6.69 

6.73 

-0.1399 

-0.1621 

-0.0699 

0.18 

4.15 

4.19 

-0.1399 

-0.1621 

-0.0699 

m 

1.44 

1.48 

-0.1399 

-0.1621 

-0.0699 

m 

-1.26 

-1.21 

-0.1399 

-0.1621 

-0.0699 

0.30 

-3.82 

-3.76 

-0.1399 

-0.1621 

-0.0699 

0.34 

-6.09 

-6.03 

-0.1399 

-0.1621 

-0.0699 

0.38 

-7.97 

-7.91 

-0.1399 

-0.1621 

-0.0699 

0.42 

-9.38 

-9.31 

-0.1399 

-0.1621 

-0.0699 

0.46 

-10.25 

-10.18 

-0.1399 

-0.1621 

-0.0699 

0.50 

-10.54 

-10.47 

-0.1399 

-0.1621 

-0.0699 

★ 


1  lb  =  4.448  newtons 
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models.  These  results  show  a  considerable  difference  between  the  unimodular  and 
bimodular  models  and  a  slight  difference  between  the  oimodular  and  multimodular 
models.  Therefore,  although  a  four-segment  model  is  a  better  approximation,  the 
two-segment  approximation  gives  nearly  the  same  results.  This  proves  that  the 
bimodular  model  precision  is  a  good  approximation. 

The  values  of  the  first  three  mode  shapes  for  the  clamped-free  case  are  pre¬ 
sented.  Finally,  the  effects  of  axial  translatory  and  rotatory  inertia  coefficients 
on  axial  force  for  a  clamped-clamped  beam  are  discussed. 

The  transfer-matrix  method  is  found  to  be  very  effective  in  terms  of  compu¬ 
tational  time  and  also  in  terms  of  the  accuracy  of  results,  which  agree  very  well 
with  the  closed-form  solution. 
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APPENDIX  A:  FITTING  MINIMIZED  CURVES  TO  THE  STRESS-STRAIN  CURVE 
1 .  Multimodular  Case 

Consider  the  nonlinear  stress-strain  curve  shown  in  Fig.  A.l.  For  any  arbi¬ 
trary  point  (et,at)  in  the  tension  region  (ej>0),  there  are  two  straight  lines 
such  that 


(c^/e^e 

[(at  -  aft)/(£t  -  eft)](e  -  e^)  + 


(A.l) 


The  equation  of  a  stress-strain  curve  is 

o(e)  =  Ke11  ;  e  >  0  (A. 2) 

where  K  and  n  are  constants  depending  on  the  material.  To  find  the  proper  "break 
point"  (£t,at),  the  area  between  the  approximated  curve  g( e)  and  the  actual  experi¬ 
mental  curve  o(e)  has  to  be  minimized.  The  mentioned  area  can  be  expressed 


, 

E 

r 

cZ  f 

A  * 

CSl  ( c)  -  o(  e)  ]de 

0 

+ 

CgoU)  -  o(£)]d£ 
h  £ 

E 

Substitution  of  equations  (A.l)  and  (A. 2)  into  equation  (A. 3)  and  performing  the 
integrations  gives 
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Fig.  A. 3  Un i modular  model . 
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A  = 


1  t  t  K  ,  t, 

o  j  £  -  ITT  (  £  ) 


n+1 


n+T 


1/  t  ,  tw  t  t\  K  r /  t \ 
2  af  J  ''£f  ‘  e  )  '  rviT^f  ^ 


n+1 


n+1 


f  -*■  i  ~  i  ~  )  ]  I  (A. 4) 

8y  searching  in  the  region  of  a  =  (0 ,etf)  x  (0 one  is  able  to  find  a  point 
(e^.a^)  such  that  A  is  minimized  locally.  Note  that  a  few  other  methods  (e.g., 
least-squares  method)  have  been  tried  but  it  turned  out  that  the  absolute  minimum 
point  was  outside  of  the  region  n. 

2.  Bimodular  Case 

For  this  case,  the  least-squares  method  has  been  used.  As  shown  in  Fig.  A. 2, 
there  is  a  line  such  that 


I  = 


:f  t  n  2 

[Efa  £  -  Ke  ]  de 


(A. 5) 


t 


can  be  minimized  in  n.  Here,  is  the  slope  of  that  line.  By  taking  the 
derivative  of  equation  (A. 5)  and  equating  it  to  zero,  one  has 


dl 
dE.  1 


*  2 


[E.  ^e  -  Ken]e  de  =  0 


By  solving  equation  (A. 6)  for  Ebt,  one  obtains 

C  t  _  3K  ,  tx0'1 

Eb  -  H+?  (ef  > 


(A.  6) 


(A. 7) 


For  example,  for  aramid-rubber  in  the  tension  region,  the  following  parameters 
are  found  (the  other  constants  are  listed  in  [10]): 

nt  ■  1.22  ;  Kt  ■  1.1  x  106  psi  ;  ef*  =  0.029 

Eb1  =  10~^  (0.029) 1  *22_1  =  0.47  x  106  psi 

An  analogous  calculation  can  be  applied  for  the  compression  side  of  the  bend, 
i.e.,  Ebc  can  be  found,  provided  that  Kc,  nc,  and  e^c  are  known. 

A  similar  approach  can  be  used  to  obtain  a  best-fit  single  straight  line 
("unimodular"  approximation)  [5]j  see  Fig.  A. 3. 
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APPENDIX  8:  THE  BEAM  STIFFNESSES  FOR  RECTANGULAR-SECTION  BEAMS  OF  MULTIMODULAR 
MATERIALS 

For  the  assumed  four-segment  model,  there  are  two  different  bending  cases  in 
general,  convex  downward  and  concave  downward  bending.  In  convex  downward  bending, 
the  top  layer  of  a  beam  is  in  compression  and  the  bottom  layer  in  tension. 

A  stress  distribution  for  convex  downward  bending  is  shown  in  Fig.  2.  As 
one  might  notice,  the  location  of  zp,  ac  and  at  fall  within  the  depth  of  the  beam 
(this  is  case  is  the  most  general  case).  Substitution  of  equation  (1)  into  equa¬ 
tion  (8)  and  using  equations  (3),  (4),  and  (5)  leads  to 


rac 


N  = 


[<E1  (ac-zn)  +  <E2  (z-ac)]dz  +  j  <E]  (z-zp)dz  + 
-h/2  a„ 


z_  at 

n  r  t 


<E1t(z  -  zn)dz 


h/2 


[<E1t(at-  zn)  +  <E2t(z  -  at)]dz 


(B.l) 


and 


rZn 


-h/2 


[<E1c'(ac- zn)  +  <E2  (z-ac)]z  dz  +  <E]  (z-zp)z  dz 


a 


h/2 


+  <E1t(z-z|l)z  dz+  [<E1t(at  -  zp)  +  <E2fc(z  -  at)]z  dz 
zn  at 

Equations  (B.l)  and  (B.2)  can  be  written  in  the  following  form 

a*  h/2 


N  =  (-<zn){[ 


E2cdz  + 


rZn 


-h/2 


+  C- 


E-jCdz  + 


’E^dz  + 


E2tdz3 


c  n 

a. 


h/2 


h/2 


E2cdz  + 


-h/2 


E-jCdz  + 


-h/2 


E^dz  - 


E^dz]} 


(B.2) 
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fac 

+  (<)([(  E2cz  dz 


2  a  h/2 

'  n  c  r  t  i  *. 

+  I  E]  z  dz  +  ;  E|Z  dz  +  j  Eg  z  dz 


,ac  ,ac  h/2  h/ 2 

+  [-l  E2cacdz  +  j  E1cacdz  +  |  E^dz-f  E^dz]: 
-h/2  -h/2  a.  a+ 


(B.  3) 


a„  z„  a*  h/2 

r  c  r  n  r  rt*.  /■  * 

M  =  (-<zn){[  E2cz  dz  +  j  E-j cz  dz  +  j  E^z  dz  +  |  EgZ  dz] 


ac  ac  h/2  h/2 

+  [-  E2cz  dz  +  E^cz  dz  +  j  E^z  dz  -  j  E^z  dz] 
-h/2  -h/2  at  & 

a  z  a.  h/2 

+  (<){[]  E2cz2dz+[  E1cz2dz  +  |  E1tz2dz+f  E2tz2dz] 


c  ac  h/2  h/2 

+  [-j  E2cacdz  +  |  E1cacdz+|  E^dz-f  E^dz]} 


(B.4) 


Combining  equations  (7)  and  (9),  one  gets 
N  =  (-<z  )A'  +  <B ' 


M  =  (-kz  )B"  +  <D‘ 


Comparison  of  equations  (B.3)  and  (8.5)  with  equations  (B.4)  and  (B.6)  and 
considering  equations  (10),  one  finds  that 


(B.5) 

(B.6) 


:A  =  (Ej  -  E2/dz  +  J  (E-|  -  E2  )dz 
-h/2  af 


(E,0  -  E2c)acdz  +  f  (E^-E^Ja^z 
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ac  h/2 

CBM  =  j  (E^-  E2c)z  dz  +  |  (E1t-E2t)zdz 
-h/2  kc 


h/2 


(E1c-E2C)acz  dz  + 


(E,  -  E, 


')at2 


dz 


(B.7) 


-h/2  at 

Eight  possible  cases  may  occur  depending  on  the  location  of  zn>  ac, 
and  a^.  These  cases  have  been  analyzed  as  the  same  as  the  general  case  as  follows 
(for  convex  downward  bending). 

There  are  seven  more  possible  cases  which  have  been  discussed  in  Ref.  [5]. 
However,  the  case  explained  here  is  the  most  general  case. 


APPENDIX  C:  TRANSFER-MATRIX  FORMULATION 


Under  harmonic  excitation,  the  steady-state-response  displacements  u,  w,  and 
„  are  also  harmonic  in  time.  Therefore,  equations  (11)  can  be  written  as 


N  =  -  n2Pu  -  «2R(p  ;  Q  =  -  n2Pw  -  q(x,t) 

>  A  ,  X 

M  -  0  =  -  n2Ru  -  n2 lip 
>  X 


(C.l) 


where  all  of  the  barred  quantities  are  amplitudes,  i.e.,  N(x,t)  =  N(x)  sin  nt,  etc. 
The  continuity  at  each  station  implies 


=,R  ■  =,L 


”iR  =  "iL 


-  R  _  -  L 

*i  ‘  *i 


(C.  2) 


(R  and  L  denote  right  and  left,  respectively) 

Also,  equation  (C.l)  in  finite-differential  form  for  each  station  i  is 
=  NiL  -  n2PG.L  ;  Q.R  =  q.L  -  n2Pw.L  -  Q. 

=  M.L  -  o2RuiL  -  n2I^.L 

where  is  the  concentrated  load  amplitude  at  station  i. 

Equations  (A. 2)  and  (A. 3)  are  written  in  matrix  notation  as 


(C.3) 
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(C.4) 


or 

R  L 

CS]  =  [T_]  [S] 
i  1  i 


where  [T  ]  represents  station  matrix  at  station  i.  In  matrix  notation  the 
s  i 

equilibrium  equation  for  each  field  under  a  distributed  load  q(x)  is 
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where 

y  =  B'B"  -  A'O' 


At 

qU)dc 

o 


At 

5qU)dc 

0 


(C.6) 


Values  of  Km  and  Kq  for  various  loadings  are  listed  in  Table  5.  Equation  (C.6) 
also  can  be  written  as 
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30 


L  R 

[S]  =  [T,]  [S] 


i+1 


( C .  7) 


l  1 


The  matrix  [T . ] .  is  called  the  field  matrix. 

J  * 

APPENDIX  D:  COMPUTATION  OF  zp,  ac»  AND  at 

For  multimodular  beams,  the  following  equation  is  not  sufficient  to  deter¬ 
mine  the  neutral -surface  location  z_ 


_  B'M  -  D'N 
zn  “  A'M  -  B"N 


(D.l) 


even  for  cases  where  N  =  0 

z  n  =  B*/A'  (D.2) 

Two  more  equations  are  needed  for  computing  zn  because  the  stiffnesses  are  not 
only  dependent  on  zp  but  they  are  functions  of  ac  and  at  as  well.  By  comparison 
of  Figures  1  and  2,  one  can  get  (for  the  convex  downward  case) 


ac  =  (£lC/£fC)  (h/2  +  zn)  -  zn  (D.3) 

at  =  (eiVe/)  (h/2  -  zn)  +  zp  (D.4) 

For  the  concave  downward  case 

ac  =  (€.,C/cfC)  (h/2  -  zn)  +  zn  (D. 5) 

at  =  (h/2  +  2n)  ‘  zn  {D’6) 


The  system  of  nonlinear  equations  (D.l),  (D.3),  and  (D.4)  for  the  convex  down¬ 
ward  case,  or  equations  (D.l),  (D.5),  and  (D.6)  for  the  concave  downward  case, 
can  be  solved  by  using  iteration  of  the  Gauss-Seidel  type. 
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